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Abstract 

With a bottom-up approach, we consider naturalness in the MSSM and NMSSM. Assuming the 
light stops, the LHC gluino search implies that the degree of fine tuning in both models is less than 
2.5%. Taking the LHC hints for the SM-like Higgs boson mass rrih ~ 125 GeV seriously, we find 
that naturalness will favor the NMSSM. We study the Higgs boson mass for several scenarios in 
the NMSSM: (1) A large A and the doublet-singlet Higgs boson mixing effect pushing upward or 
Q-I pulling downward m^. The former case can readily give the di-photon excess of the Higgs boson 

decay whereas the latter case can not. However, we point out that the former case has a new large 
fine-tuning related to strong A— RGE running effect and vacuum stability. (2) A small A and the 
mixing effect pushing m^ upward. Naturalness status becomes worse and no significant di-photon 
excess can be obtained. In these scenarios, the lightest supersymmetric particle (LSP) as a dark 
matter candidate is strongly disfavored by the XENON100 experiment. Even if the LSP can be 
a viable dark matter candidate, there does exist fine-tuning. The above naturalness evaluation is 
based on a high mediation scale for supersymmetry breaking, whereas for a low mediation scale, 
fine-tuning can be improved by about one order. 
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I. INTRODUCTION AND MOTIVATIONS 



Supersymmetry (SUSY) solves the gauge hiearchy problem in the Standard Model (SM) 
naturally. In the supersymmetric SMs with i?-parity, gauge coupling unification can be 
achieved, which strongly indicates the Grand Unified Theory (GUT). Also, a surprising gift 
is that the lightest supersymmetric particle (LSP) services as a cold dark matter (DM) 
candidate. The recent LHC experiments and DM search experiments lead to the following 
implications: 

• The CMS [1] and ATALS [2] SUSY searches constrain the colored particles stringently 
and push their masses towards the TeV region, except that the light stops ~ 0(200) 
GeV can still be allowed [3, 4]. Light stops are good for naturalness by virtue of the 
absence of enormous fine-tuning from the stop radiative corrections on the Higgs soft 
mass square 

2 3/^(Tr4+m^) M a 
m H u £-2 log m ¥ , (1) 

where M is the SUSY-breaking mediation scale. Moreover, they maintain the discovery 
potential of colored supersymmetric particles (sparticles) at the early LHC run. 

• In the Minimal Supersymmetric SM (MSSM), the lightest CP-even Higgs boson mass 
is less than the Z-boson mass mz, and is lifted by the top and stop loop corrections. 
However, the recent Higgs search results at the ATLAS and CMS experiments imply a 
SM-like Higgs boson with mass around ~ 125 GeV [5] , which inspired extensive studies 
within the MSSM [6-8]. Such a relatively heavy Higgs boson typically requires TeV 
scale stops and an anomalously large \A t \. In the gauge mediated symmetry breaking 
(GMSB) scenario, a large \A t \ is viable only when the messenger scale is very high [7], 
which is a prediction in the "asymmetric" gravitino scenario [9]. 

• The discovery of the relatively heavy Higgs boson and the null XENON100 experimen- 
tal result [10] place the neutralino LSP DM candidate in an unsatisfying position [8]. 
Now the LSP abundance tends to be far above the WMAP-measured value except for 
very fine-tuned coannhilations or for a Higgsino-like LSP which is, however, strongly 
disfavored by the XENON100. Note that it can constrain on the /x parameter and the 
naturalness thereof. 

In summary, the naturalness is quite pessimistic in the mSUGRA-like MSSM, although it 
is still far from the failure allegation. 

Taking all these hints, what is the most natural SUSY model one can find? The next- 
to-MSSM (NMSSM) stands in the foreground. Historically, it was proposed as a simple 
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solution to the fi— problem, which is also a naturalness problem in the MSSM (For details, 
see Ref. [11] and references therein.). As a by product, it has an advantage in increasing 
the tree-level Higgs boson mass. Therefore, it allows the lighter stops and weakens the 
correlation between mz and m^. This property not only resolves the LEP crisis, but also is 
likely to be the most promising savior of natural SUSY at the LHC (For recent discussions 
on the NMSSM with a relatively heavy Higgs boson, see Ref. [12-15].). 

In this work, we study the naturalness via mz or m^ in the MSSM and NMSSM using 
a bottom-up approach. Given the light stops and gluino with respectively the LHC lower 
bounds ~ 200 GeV and 600 GeV [3], the least degree of fine tuning involving m z is due to 
the gluino, and roughly 2.5% for both models. However, taking the hints for rrih ~ 125 GeV 
seriously, we show that the NMSSM is more natural. We study several scenarios which may 
give a relatively heavy SM-like Higgs boson in this model: 

• A large A and the doublet-singlet Higgs mixing effect pushing m^ upward, which has 
a heavier Higgs boson and an significant di-photon excess. However, this scenario has 
a new large fine-tuning due to vacuum stability and large A— renormalization group 
equation (RGE) running effects. 

• A large A but the mixing effect pulling m^ downward, A ~ 0.7 and tan/5 ~ 2 allows 
rri/ L ~ 125 GeV. An essential difference between this scenario and the above previous 
is the absence of a lighter Higgs. In addition, it is very difficult to give the significant 
di-photon excess. 

• A small A and the mixing effect pushing m^ upward. We do not have significant 
di-photon excess, and the naturalness status becomes worse. 

In all the above scenarios the neutralino LSP DM candidate is strongly disfavoured by the 
XENON100 experiment. Even if the neutralino LSP DM is fine, there still exists fine-tuning. 
Note that the above analyses are based on the mSUGRA-like model with the mediation scale 
M = Mqut; but when M is sufficiently low the naturalness can be improved by about one 
order. 

This paper is organized as follows. In Section II, we make a detailed analysis on the 
naturalness implication via mz and m^ in the MSSM and NMSSM. We discuss the re- 
lated phenomenological consequences in Section III. The Section IV is the discussion and 
conclusion. Some necessary and complementary details are given in Appendices A and B. 
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II. THE ROAD TO THE MOST NATURAL MSSM AND NMSSM 



The LHC is testing supersymmetric models. The most predictive model such as the 
Constrained MSSM (CMSSM) has been pushed to the multi-TeV region directly by the 
SUSY search or indirectly by the Higgs search. Naturalness is seriously challenged there, 
nevertheless the NMSSM still can be natural. In this Section, we will study the origin of 
fine-tuning in the MSSM and NMSSM via the bottom-up approach. 



A. Light stops and gluinos: The natural soft SUSY spectrum for rnz 

We will first discuss the naturalness implication on the Higgs sector with successful elec- 
troweak (EW) symmetry breaking defined at the weak scale, and then use RGEs to extrap- 
olate relevant soft parameters to the UV boundary and examine the naturalness in terms of 
fundamental soft parameters. 



1. Naturalness of the electroweak Higgs sector 



The Higgs parameters not only determine the EW-scale ~ rnz but also have close relation 
with the SM-like Higgs boson mass rrth- In the MSSM, the tree-level Higgs quartic coupling 
is uniquely determined by the SU{2) L x U{1)y gauge couplings via D-terms, which implies 
the tree-level upper bound on the lightest CP-even Higgs boson mass: m 2 < m| cos 2 2/3. 
To make up the gap between rrtz cos (3 and the hinted mass ~ 125 GeV, a substantial 
radiative correction from the top-stop sector is necessary, i.e., 

(2) 

with X t = A t — fi cot (3 the stop mixing and mj the geometric mean of two stop masses. 
To sufficiently lift m^, a rather heavy stop sector with m-j- > 700 GeV is needed, even in 
the maximal mixing scenario [12]. In turn we are forced to accept a large —fn 2 H u which is 
implied by Eq. (1). 

Therefore, in the MSSM such a heavy Higgs boson is obtained at the price of a fine-tuned 
Higgs sector. This can be seen from the following tadpole equation which determines the 
weak scale rriz 




in 



V , (3) 



2 tan 2 - 1 

where the parameters are defined at the EW scale. If all the Higgs parameters lied around 
m z (tj 2 ? is an exception whose contribution to the Z boson mass may be suppressed by 
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large tan 2 (3 [41]), then the determination of mz would be natural. Otherwise, one can adopt 
the following quantity to measure the degree of fine-tuning associated with mz [17] 



A z = maxjFi, 



dlnmz 



d In pi 



(4) 



with pi the fundamental parameters. The above definition can be applied to any quantity 
which is a consequence of cancellation. The putative m^ ~125 GeV renders Az < 0.1% or 
even worse [12]. 

Turn our attention to the NMSSM, in which the impact of a heavier on the naturalness 
of mz can be abated considerably, since the NMSSM specific effects are capable of lifting 
mh without heavy stops. Nevertheless, /i as well as other Higgs parameters themselves may 
hide new fine-tunings, when we are committing ourself to find a relatively heavy m^. To 
investigate the actual naturalness status of the NMSSM (here only the Z3— NMSSM is under 
consideration), we start from the Higgs sector 

W = XSH,-H d+ ^, 
-C„„ =m 2 „JH u \ 2 + ml d \H d \ 2 + m||S| 2 + [\A X SH U ■ H d + ^A K S 3 + fc.c.) . (5) 

On a proper vacuum we have v s = (S) ^ 0, therefore the li = Xv s is generated dynamically. 
The non-observation of charginos at the LEP2 gives a bound \i > 104.5 GeV. To demonstrate 
how new fine-tuning may arise in determining /i, analogously to Eq. (3) we trade the order 
parameter v s with \x in the singlet scalar tadpole equation, and get 

/ K 2 K \ 1 

ix 2— /i 2 + -A K y, + m 2 s + A 2 ^ 2 - Xkv 2 sin 2/3 - -X 2 A x v 2 sin 2/3 = 0. (6) 
V A A J 2 

Then if all the NMSSM specific soft parameters lie much above the weak scale, a small \x 
will be a result of fine-tuning. 

A remark deserves attention. Tied to \x and hence the weak scale directly, the chargino 
mass is a key to understand the naturalness of SUSY. The absence or the discovery of a 
light chargino ~ 200 GeV in the future will shed light on it. Interestingly, we will find that 
the naturally heavier m^ and smaller ii are inherently consistent in the NMSSM . Concretely 
speaking, \x will be found to automatically fall into the narrow region 100-300 GeV, which 
ensures the proper doublet-singlet (H/ S) mixing as well as a non-tachyonic light stop secctor. 

A subtle hidden fine-tuning is associated with vacuum stability. It arises in the large A 
limit, when the strong RGE effects significantly change the values of soft parameters, saying 
m|, during the RGE flowing. We postpone to Section II B 3 for the concrete discussions on 
how does it render the fine-tuning in vacuum stability. 
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Finally, we have to remind that the (3— angle is dynamically determined through the third 
tadpole equation 

Sm 1 m\ u + m\ d + 2/i 2 + X 2 v 2 ' 1 ' 

which seemingly does not invoke fine-tuning. However, at least in the MSSM this equation 
can be transformed to the determination of the CP-odd Higgs boson A mass square: 

m\ = 2Bfi/ sin 2/3 = m 2 Hu + m 2 Hd + 2/i 2 + A V. (8) 

If rriA were measured to be far below the weak scale, then it would introduce a new source 
of fine-tuning comparable to the one of mz ■ 



2. Feed back to the boundary 

In the previous subsection, we only outline a natural Higgs sector at the weak scale. But 
its parameters are not fundamental and should be evolved to the mediation scale M. In 
this work two prevalent mediation scales are considered: the GUT scale Mqut mediation 
with M = Mqut in the mSUGRA-like model, and the low mediation scale with M = 10 5 
GeV in the gauge mediated SUSY breaking (GMSB) (For a review, see Ref [19].). Since 
we are endeavoring to explore a less fined-tuned model, some unification assumptions in the 
mSUGRA model are dropped. Concretely speaking, the Higgs sector parameters are free, 
and the stops are treated in a special role just as in the SUSY framework with the light 
third family sfermions [4]. Howbeit, the gaugino mass unification is preserved, and similar 
strategies will be adopted in the GMSB. 

SUSY with high mediation scale suffers large fine-tuning, and one of the main object 
of this work is to explore how natural it can be. The exploration relies on one fact: the 
soft parameters at the low energy SUSY-scale Ms = 1 TeV are polynomial functions of the 
initial values at the GUT scale, with coefficients being functions of the Yukawa and gauge 
couplings [20]. For example, the Higgs doublet soft mass squares receive the large corrections 
from gluinos and/or squarks: 

m 2 Hu ~ - MA 2 + A9A t M 1/2 - 2.l7M 2 /2 - .llm 2 Hd + A2m% u - .08m| - .33m| c - .43m|, (9) 
m 2 Hd -0.02A 2 - mA t M 1/2 + .50Mf + .73m^ d - .14m^ - .13m|, (10) 

where the parameters with bar are defined at the UV boundary. We have set A = 0.62, k = 
0.30, tan (3 ~ 2.5, m t ~ 173.4 GeV at Ms (it is in the bulk region studied numerically later), 
and these inputs will be adopted throughtout this Section. 
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Now we are able to evaluate the fine-tuning of mz with respect of the various fundamental 
soft term parameters. The Z— boson mass square is expressed as 

2 

^ ~2.68M 1 2 /2 - .53rh 2 Hu + .27m 2 Hd 

+ .53m| + .40m| c - .06A 2 - .22A t M 3 + .07m|. (11) 

It is supposed that fi ~ 100 GeV does not contribute much to the fine-tuning of mz, so it is 
not included here. Then the possible large tuning, in light of Eq. (4), are calculated to be 

F Ml/2 ~ 2 x 2.68Mf/m|, = 1 x 0.53m 2 H Jm 2 z , F ij = 2x 0.06A t 2 /m|, 

= 1 x 0.40m| c /m|, = 1 x 0.53m|/m|, F m | = 1 x 0.06m|/m|, ... (12) 

where the dots denote irrelevant terms. Each F, the degree of fine-tuning, is the fine-tuning 
coefficient (e.g., 0.53 for F m 2 ) which encodes the RGE effect, times the initial value over m|. 
Note that we are regarding m\ rather than rfiH u as fundamental. Thus, the corresponding 
tuning is reduced by half. Otherwise, the value should be doubled. 

Now we investigate the naturalness implications on various SUSY breaking soft parame- 
ters. The gluino mass tends to be the dominating source of fine-tuning, as a consequence of 
the long RGE running from the GUT-boundary . The LHC SUSY search places a bound on 
the gluino mass Mg > 600 GeV corresponding to the initial value iWi/ 2 « 250 GeV (fixed 
hereafter), which gives the least degree of fine tuning: 

A z > F Ml/2 ~ 40. (13) 

Throughout this work we will take it as a referred value of the degree of fine-tuning for the 
natural NMSSM, unless otherwise specified. Note that from Eq. (7) we obtain 0.27m 2 Hu + 
0.62m^ d ~ 2-B/i/sin2/3 ~ 10 5 GeV 2 , where B ~ A\ ~ 400 GeV. So typically we have 
fn 2 Hu > ^H d (~ "frig) ~ 10 6 GeV 2 , which implies that the m^— term usually contributes the 
largest tuning to get the weak scale. 

In the light stop framework, it is supposed that stops do not lead to large fine-tunings. 
However, the multi-TeVs A t also results in an enormous fine-tuning and roughly gives the 
bound | At | < 1.7 TeV to satisfy F^ t < F Ml/2 . On the other hand, it is likely that a large 
| At | is necessary to achieve a large mixing in the stop sector, which is important to enhance 
rrifi with stops as light as possible. To see that clearly, we consider the following equations: 

A t ~0.19A t - 0.06A A - 1.83M 1/2 , (14) 
m| c ~0.61m| c - 0.26m| - 0.04A 2 + 0.14A t M 1/2 + 3.04M 2 /2 - 0.2m% u + 0.03m|, (15) 
m| ~0.84m| - 0.13m| c - 0.02A 2 + 0.07A t Mi /2 + 4.33M 2 /2 - 0.13m^ + 0.02m|. (16) 
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As one can see, — A t ~ 0(1) TeV helps to achieve a large stop mixing by not only giving a 
large A t but also reducing the stop mass squares. 

Comments are in orders. In the first, it can be explicitly seen from Eq. (3) that the m 2 z 
strong gluino mass dependence from m 2 Hu (as well as the weaker dependence from mj^J is 
exacerbated as tan /3 becomes smaller, so a small tan (3 < 2 is not favored for getting the 
optimal rriz- Next, the previously discussed fine-tunings involving rriz originate from the 
stop and gluino sector, but are independent on the singlet sector. So the conclusions can be 
applied to both the MSSM and NMSSM. 

Due to the large A, the soft terms in the NMSSM singlet sector develop strong dependence 
on the MSSM terms, e.g., the Higgs soft terms. That may give rise to new fine-tunings. The 
singlet soft terms are given by 

m 2 s ~ - 0.30rh 2 Hd - 0.21fn 2 Hu + 0.10m| + 0.25m| + (17) 
A\ ~ — 0.24A t + 0.31A A + 0.27M 1/2 , 

A K «0.3L? t + 0.814, - 1.06 A x + 0.06M 1/2 . (18) 

Several remarks about the RGE effects on the low energy soft terms are in orders: (A) \ms\ 
develops a strong dependence on \mu ud \ ~ 1 TeV, and thus it is typically close to the TeV 
scale in the absence of fine-tunings. On the other hand, in Section II B 3 we will find that 
vacuum stability typically requires \m s \ < 100 GeV. As a consequence of this tension, m| 
usually renders a rather large fine-tuning 

A s ~ 0.25m|/m|. (19) 

(B) The fine-tuning coefficient of m| is self- decreased by a large k. Varying k from 0.1 to 
0.3, it decreases about three times (with A fixed). We will use a simple function to fit this 
behavior, which is good enough to calculate As numerically. (C) Similarly, the A\ and 4 
components in A K and 4 are sensitive to A and k. Especially, in the limit A — > 0.7, the Ax 
component increases much in A K whereas self-squeezes substantially in A\. Therefore, the 
expected natural order is \A K \ > \A\\, otherwise a new fine-tuning may occur. 

Now we turn the attention to the GMSB with low mediation scale. In contrast to the 
high mediation scale case, for M = 10 5 GeV the gluino effect is reduced and the fine-tuning 
is thus considerably improved. Concretely, now the Z boson mass square assumes a form 

m% ~ 2 (o.08M 3 2 + 0.17m| + 0.15m| c + 0.12A t 2 ...) . (20) 

Compared to Eq. (11), it is clearly seen that the gluino mass dependence has been reduced 
one order of magnitude, although the squark influence is still rather significant. Thus again 
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light stops with masses ~ 600 GeV, are elements for natural SUSY. As a comparison, the 
squark soft mass squares are now given by 



m| c ~0.91mJ c - 0.08m| + 0.41M 3 2 - 0.06A 2 t , 
m| ~ - 0.04m| c + 0.95m| + 0.44M 3 2 - 0.03A 
A t «0.72A t - 0.35M 3 . 



t ■ 



(21) 
(22) 
(23) 



As for the singlet sector, because of the much shorter RGE running, the A— RGE effect is 
reduced considerably or even ignorable and then not listed. 

At the messenger boundary, we have assumed that the stop mass squares are naturally 
and properly small. Thus, in principle the fine-tuning can be reduced to a level of 10% or even 
better, depending on the gluino mass (Mgr=750 GeV to get 10%). However, in the minimal 
GMSB the gluino mass is tied with the first and second family squark masses. The exclusion 
of such light squarks means that Mg should lie above the TeV scale. Consequently Az is 
not quite optimistic again. But recalling that generating acceptable large gaugino masses is 
a generic problem in the GMSB with dynamical SUSY-breaking [21], we conjecture that the 
gaugino masses may be suppressed comparing to the squark soft masses. Then again Az 
is directly related to the LHC gluino search, and Az ~ 20% is still viable from the present 
data. In a word, naturalness strongly prefers the GMSB-like models. 

To end up this section, we summarize the naturalness implications on the GUT-scale soft 
terms for high scale mediation: 

• The stops and gluino should be as light as possible. Moreover, A t enters mz and its 
rough naturalness bound is \A t \ < 1.7 TeV. 

• By virtue of a large A, one expects that the low energy soft terms show m| > O(10 5 ) 
GeV 2 and \A\\ < \A K \. Otherwise, extra fine-tunings are introduced. 

• In the weak coupling limit the A— RGE effect is reduced greatly, so the NMSSM does 
not introduce extra fine-tuning. 

If the mediation scale is sufficiently low, naturalness of both the Z boson and Higgs boson 
masses can be improved much. In light of the previous analysis, the inventory of the natural 
NMSSM in the deformed GMSB is: the properly suppressed gaugino (at least gluino) masses, 
light stops, and of course, a sufficiently low messenger scale . We leave this more optimistic 
scenario for a future work [22]. In what follows we will focus on the high mediation scale 
case. 
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B. Mixing scenario: a natural heavy Higgs in the NMSSM 



In the previous subsection we have shown that on the naturalness of mz alone, the LHC 
imposes the common lower fine-tuning bound for both models. Thus the NMSSM does not 
have an obvious advantage over the MSSM, and maybe it is even more tuned due to a large 
A. However, it is preferred when taking the Higgs boson mass into account. In the LEP era, 
the NMSSM alleviates the little hierarchy problem via two different strategies: Pushing nih 
above 114.4 GeV with a large A and a small tan /3 [23] or a proper H/ S mixing effect [24] [42], 
or alternatively allowing below the LEP bound but requiring non-standard Higgs decays 
to escape from the LEP search [27]. Now the ATLAS and CMS hints exclude the second 
scheme, but the first scheme still works and may be favored viewing from the Higgs di-photon 
excess. 

Within our knowledge, previous literatures on the naturalness of are based on the 
LEP bound, and treat the H/S mixing effect mainly using a numerical method [24, 26]. 
Thus intuitions may get lost. As one of the main object of this work, we employ a detailed 
analysis of this effect using a semi-analytical method. We start from the following two by 
two mass squared matrix 



which approximately encodes the H/ S mixing effect in light of Appendix A. We denote 
its eigenstates as H\ 2 with mass ascending. (M^ s )n gives the well-known NMSSM upper 
bound on the SM-like Higgs boson mass without taking the mixing effect into consideration. 
We are ignoring the correction from the 13— mixing effect, which decreases (M^ s ) 2 2 by an 
amount (M|)^ 3 / M% x . Since (M|) 13 oc cos 2/3. The decrease may be appreciable for a large 
tan/3, but it will not affect our main conclusion. So we leave this subtle effect to numerical 
studies. In addition, the top-stop loop correction gives 5 t — 50 GeV — 80 GeV for ~ 300 
GeV and a moderate mixing effect. Throughout our analysis, such a referred value will be 
used as the premise of defining natural SUSY. 

1. A general analysis on the H/S mixing effect 

As a general discussion on the mixing effect and its implication, we start from a general 
matrix structure rather than confining to the NMSSM. So our discussion can be applied 
to many new physics models showing a H/S mixing, e.g., the SM Higgs sector extended 
with an extra scalar singlet in the context of dark matter models [28, 29]. For an arbitrary 
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2x2 symmetric real matrix M 2 , the necessary and sufficient condition that the mixing effect 
pushes the larger diagonal element to an even larger eigenvalue (namely the SM-like 



Higgs boson mass) is 



M n >M 2 2 2 ,M 2 



M n M 2 2 > 



(25) 



The second equation is the minimum condition or vacuum stability condition. For our 
purpose, a natural Mn should be around 120 GeV. The two eigenvalues of M 2 are written 

as 



where A± = M\ x ± M| 2 . The mixing effect pushes upward the "mass" (in the flavor basis) 
of the doublet M n to the eigenvalue m 2 , whereas pulls downward the singlet "mass" to m\. 
We will use the pushing or pulling effect to describe these two faces of the mixing effect for 
short. 

The implications of the mixing effect are two-folds: One is the mixing angle between 
states, and the other is the pushing or pulling effect on "Higgs boson masses" . The unitary 
rotation diagonalizing the mass square matrix is given by 



Thus the mixing angle is determined solely by the ratio A_/2M 2 2 , which can be large even 
for a very small off diagonal element M 2 2 <C M\ x . Hence we conclude: Although the mixing 
angle is of importance in the Higgs collider phenomenology, it can not properly measure the 
amount of pushing upward or pulling downward m^. 

To estimate the size of the enhancement effect, we go to some special but well- motivated 
limits: 

Decoupling limit When the off diagonal element is so small that 2M 2 2 <C |A_|, it is ap- 
propriate to define the decoupling limit where both the mixing angle and enhancement 
is vanishing small. In this limit the doublet mass eigenvalue is given by 



Since 2M 2 2 /M 2 1 < 2M 2 2 /A_ < 0.3 by definition, the enhancement is indeed ignorably 
small in the decoupling limit: m n — Mn < l%M n . 




(26) 




(27) 




(28) 
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Marginal push Considering M\ x 3> Mf 2 , one can imagine that the mixing effect should 
be marginal for a properly large M 2 2 . In this case the heavy Higgs boson mass square 
is approximated to be 

m\ « M 2 

and the mixing angle is tan6> ~ MfjMf 2 > |Mn/M 22 | by Eq. (25). For M n = 120 
GeV, the marginal mixing effect is able to give the necessary push, e.g., M\ x = 3M| 2 
and Mf 2 = Mfj8. While M n = 115 GeV does not work. 

Significant push The large enhancement effect will happen when M\ x > M| 2 and M 2 2 3> 
|A_|. Then we have 

9 1 / 9 o \ 9 (M 2 - M? 9 ) 2 , N 

m? « - (M{\ + M 2 2 2 ) + M 2 2 + 1 u gM? 22J ■ (30) 

In principle, the mixing effect could be able to make the heavier Higgs approaching 
the maximum value ~ y/2Mn, together with a light Higgs boson. The mixing angle 
in this case is large. 

More generically, the numerical plot in Fig. 1 gives some insights into the features of H/S 
mixing, confronting with the collider bound: 

• For a fixed pushing strength (saying 1.05), the doublet fraction varies in the region 
~0. 50-0. 90. For the fraction ~ 0.5, the diagonal elements are quite degenerate, which 
implies that the two eigenvalues are roughly equal. In addition, such a large mixing 
renders the definition of a SM-like Higgs boson not very clear. And generically the 
h — > 27 rate is reduced and then not favorable. If the fraction is 0.9, the heavier Higgs 
field h = if 2 pairs with a lighter singlet-like Higgs boson Hi, whose mass is typically 
below the naive LEP upper bound. But it is likely that at some intermediate value of 
the doublet fraction, aided by a large enough M 22 , the lighter Higgs Hi is allowed by 
the LEP. 

• We can also consider the pulling scenario in which the large A and small tan /3 have 
already lifted mh up to ~ 125 GeV. Thus, we have to preclude the pulling effect from 
decreasing too much. This scenario is phenomenologically distinctive owing to the 
absence of a new light Higgs boson other than h = Hi. 

In the following sections, we will explore the natural Higgs scenarios along these lines. 
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FIG. 1: The contour plot of the enhancement factor m\jM\\ (thick black lines, m\jM\\ = 
1.05, 1.10, 1.15, 1.20 respectively from bottom up) and sin 2 9 the doublet component of h (dashed 
lines) in the M22/M11 — M12/M11 plane. 

2. Large A: a no-go for the decoupling limit 

The NMSSM with a large A ~ 0.7 and a small tan/3 can lift to 125 GeV even 
without a mixing effect and large loop-correction. However, it does not guarantee a natural 
NMSSM. As argued in the previous subsection, tan/3 ~ 1.7 mildly exacerbates the fine- 
tuning involving rrtz- On top of that, a new fine-tuning may emerge when we take the RGE 
effects into account. 

A large A tends to pull the SM-like Higgs boson mass rather than push. To show this, 
we rewrite M\ 2 which is given in Eq. (24) as 



If there is no substantial cancellation, the first term 2A/xt> alone shows that M| 2 > M\ x is 
needed to avoid tachyon even if /x tracks the lower bound ~100 GeV. Consequently Eq. (25) 
is violated and the mixing effect pulls down the SM-like Higgs boson mass. 

To minimize the pulling effect, we resort to the decoupling limit. It requires a quite 
heavy singlet sector, which will be found to be impossible. To see it, we consider the explicit 




(31) 
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expression for Mf 2 : 

M 2 2 2 = A sin 2/3 + 4^ 2 + ^u. (32) 

At first, the large A and putative small /i negate the possibility of the second term to exceed 
the weak scale much. Next, if we count on A\ sin2/3/2/i 3> 1, from Eq. (31) and Eq. (36) 
we know that it will render the CP-even Higgs boson mass matrix tachyonic. Finally, the 
CP-even and CP-odd singlet mass squares receive opposite contributions from the kA k term 
(see Appendix A), so the kA k dominance in M| 2 is also excluded. Eventually we affirm 
the no-go in the natural NMSSM with a large A and whatever value of tan /3: A heavy and 
decoupling singlet sector is impossible and large kA k is excluded as well. In other words, a 
proper cancellation is inevitable, and we will discuss the naturalness of this cancellation in 
the following. 



3. Large A: push versus pull 

Previously, it is found that although a large A effectively enhances the Higgs boson mass, 
cancellation is necessary. So we have to contemplate whether or not the cancellation is 
tolerated by naturalness. In light of Eq. (31) we define 

C A = f^ + f)sin2/?, (33) 



2/i A, 

which should be quite close to 1 to sufficiently reduce M\ 2 through cancellation. To maintain 
perturbativity up to the GUT-scale, A > 0.6 means that k should be moderately smaller 
than A. Then A\/2fi plays the primary role for cancellation (But k also plays an important 
role if this ratio is relatively small). So we have 

Ax - (34) 
sm2p 

As one can see, the absence of light charginos requires a large A\, especially when tan /3 is 
relatively large. 

The parameter space is quite predictive once naturalness and Higgs boson mass conditions 
are imposed. First of all, tan f3 ~ 2 and A ~ 0.7 are necessary for a relatively heavy SM-like 
Higgs boson. Next, we require /i ~ 200 GeV for the sake of significant mixing effects and 
naturalness. Then A\ falls into a narrow region around 400 GeV due to Eq. (34). Those 
observations are confirmed by Fig. 2. Finally, n should be properly large so as to suppress 
or forbid the invisible decay of Higgs into a light neutralino pair. 
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The aforementioned cancellation itself does not mean fine-tuning, since Ca can be deter- 
mined dynamically by the singlet tadpole equation Eq. (6). From this equation we get 

mf + 2(k/A) V + (k/X)A k h 
C A -1 = ^ , (35) 

where Xv ~ 100 GeV. Therefore, it seems that the right-handed side can be naturally at 
~ 0.1 due to a sufficiently small \ms\, with the other terms suppressed by small enough k 
and A K . But this does not address the naturalness problem correctly if the RGE effects are 
taken into account. In fact, \mg\ < 30 GeV leads to fine-tuning worse than 1%, which is 
manifested in Eq. (17). So the really natural \mg\ should be properly large, e.g., \mg\ ~ 10 2 
GeV. As a result, Eq. (35) implies a larger 2(k/A) 2 /U 2 , which forces k ~ A (So the PQ-limit 
is disfavored). On the other hand, a larger k further improves naturalness by squeezing the 
fine-tuning coefficient of frig. These points will be confirmed numerically. 

Another source of fine-tuning may creep in. We have shown that a large \A\\ and a small 
\A K \ is favored to avoid tachyonic Higgs states. Thus, we expect the order \A\\ > \A K \ for 
vacuum stability. However, in Section II A 2 it was found that the natural order should be 
reversed in the large A scenario. Consequently it gives rise to a source of tuning. Fortunately, 
even taking \A t \ < 1.7 TeV into account, naturalness merely places a loosely lower bound 
on \A K \. For example, setting A x = 400(600) GeV and A t = -1.15 TeV, then A Ak = 40 
means \A K \ < 30 (70) GeV, which is always satisfied in our numerical study. 

The mixing effect on rrih depends on the singlet mass square M| 2 , whose different values 
lead to different scenarios of lifting the Higgs boson mass. We rewrite M| 2 using Eq. (35) 



M 2 2 



,•2 



* (C A - j sin 2/3) A V + 4^/x 2 + jA K ». (36) 

There may be a further small correction due to the 13-mixing effect which reduces Mf 2 . The 
first term is positive and < M^. The other two terms also play some roles, since k can not 
be very small. Then, two different cases arise: 

• The pushing region is realized when h = H 2 and the lighter Hi is singlet-like. The 
LEP bound on mu 1 means Mf 2 should be properly large, which can be seen in the 
top-left plot of Fig 2. In the small A\ region where the first term of M| 2 is small, a 
large k is required. From Fig. 2 we see that \i is automatically bounded by fi < 240 
GeV (The following pulling scenario has a similar property). There are two reasons to 
understand the smallness of fi, both related to the heavy Higgs: (1) A small fi is good 
for sufficient mixing; (2) We need light stops, a large mixing — A t as well as a smaller 
tan {3, and then X t would be so large that the color symmetry breaks if fi ^> mz- 

• k/\ is relatively large (~ 0.4, see Fig. 2) and the negative third term does not cancel 
the first two terms. Then M| 2 > M\ x and we get the pulling scenario. There, asides 
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from the relatively heavy SM-like Higgs boson h = Hi, exists a mildly heavier singlet- 
like Higgs boson H 2 . We do not need to worry about the LEP bound, but should keep 
M\ 2 sufficiently small to minimize the pulling effect. This can be seen by comparing 
the distribution of Ca in Fig. 2, where the pulling scenario has much more points 
crowding around Ca = 1 than the pushing scenario. Note that the pulling scenario 
usually suffers more fine-tuning with a reason traced back to Eq. (35), where a rather 
large 2(k/A) 2 / u 2 has to be canceled by a larger \m 2 s \ and a larger m 2 H thereof. 




140 160 180 200 220 240 140 160 180 200 220 240 260 

H(GeV) H(GeV) 



FIG. 2: Predicative natural NMSSM parameter space. Top-left: tan j3 — A\ (red) and tan /3 — pt 
(green) plots in the pushing scenario. Bottom-left: Plots of the three ratios as functions of /U in the 
pushing scenario. Top-right and bottom-right: Corresponding plots in the pull scenario. A = 0.60 
is fixed, while the other parameters and physical constraints are given in Section III A. 



4- The small A scenario 



As mentioned in subsection II B, a mixing effect only may be able to push the SM-like 
Higss boson mass close to 125 GeV. This happens in a distinctive parameter space where A 
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is moderately small and then the tree-level Higgs boson mass reduces to the MSSM case, so 
a large tan (3 is required. 

To get an overall impression on the feature of the parameter space of this scenario, we 
again follow Eq. (25). First, M 22 is approximated to be 

M 2 2 2 w4(k/A) V + (k/X)A kI i. (37) 

In M 22 a term proportional to (\v/2fi) M\ 2 has been ignored, on the ground of M\ 2 < M\ x 
and a small A. M 22 should be around 100 GeV for a significant push. Recalling that a large 
kA k is inconsistent with vacuum stability, the unique option for M22 is 

/x n /100GeV\ 
k/\ ~ 0.5 ( J . (38) 

This constraint on k has important implication on the singlino mass and then the LSP phe- 
nomenology. Now the non-diagonal element Mf 2 = 2A/xw(l — Ca) can be readily suppressed 
due to a small A, even without turning to a small 1 — Ca via cancellation. In addition, again, 
like the pushing scenario with a large A, here a small [i is required to control the sizes of 
Mf 2 and M 2 2 2 . 

Although the small A scenario requires rather large cancellations between the Higgs pa- 
rameters to satisfy the tadpole equations which contain X 2 v 2 , they are not new sources of 
fine-tuning. As mentioned in Section II A 1, such cancellations are dynamically implemented 
and do not determine any particle mass scale. Actually, in the small A case the fine-tuning 
associated with m| disappears, since the strong A— RGE effects on m| are removed. 

III. THE NUMERICAL STUDIES IN THE NATURAL NMSSM 

The natural NMSSM provides a quite attractive framework from the point view of discov- 
ery potential, because of the lightness of neutralinos, the third family squarks, gluino, as well 
as Higgs boson. In this Section we present numerical studies with the aid of the NMSSM- 
Tools 3.2.0 [30], not only affirming the points developed previously but also providing some 
important phenomenological consequences. 

A. The Higgs boson decay: naturalness and 27 excess 

First, we study the large A scenario. Before the numerical study, we present the setup 
for the parameter space in the natural NMSSM. We fix the stop sector at Ms following the 
criterion of minimizing the fine-tuning from mz'- 

m ~ = 350 GeV, m Dc = 300 GeV, A t = -500 GeV, (39) 
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which leads to the lighter stop mass mr ~ 200 GeV. In terms of Ref. [3], the current 
LHC data still allow such a light stop, and we will take this lower bound on the stop sector 
hereafter. Of course, one can choose other configurations, but our choice is likely to approach 
the least fine-tuned stop sector with a realistic spectrum. Other sfermions are assumed to 
be irrelevant. Although this assumption does not hold water for extremely heavy sfermions 
which significantly modify RGEs, we suppose that such effects are sub-leading. To illustrate 
the sensitive dependence of the degree of fine-tuning on A, two cases with A = 0.60 and 
A = 0.62 will be studied for comparison, whereas k, tan/3, /i, A K and A\ vary freely. We 
keep the points satisfying all constraints given in the NMSSMTools except the WMAP and 
XENON 100 bounds which will be studied specifically. In addition, we require the SM-like 
Higgs boson mass rrih > 120 GeV and its doublet fraction should be larger than 0.8. 
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FIG. 3: The degree of fine-tuning As and F m 2 as functions of k (top) and m/, (bottom). Left 
panels: A = 0.60; Right panels: A = 0.62. Naturalness favors larger k. A heavy Higgs tends to be 
at the price of larger fine-tuning. 

To calculate the degree of fine-tuning, we record the crucial soft parameters at the GUT- 
scale boundary, which are obtained by RGE evolving the parameters from Ms to Mqut- 
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The calculation is based on some observations and reasonable approximations: (1) F m 2 
is the largest degree of fine-tuning involving mz (In Eq. (11) m\ plays the major role to 
cancel My 2 ), if it exceeds Fm 1/2 - We have the approximation 



0.42 x tan 2 (3 + 0.15 m l H . 
F-2 w — — f40l 

± mi, -, , 9 a 2 ' V^"/ 

ff « 1 — tan 2 p m z 

which varies only slightly for A = 0.62 and A = 0.60. (2) As discussed in Section II A 2, As 
defined in Eq. (19) is much more complicated than F m 2 , and we use a function of k to fit 

H u 

the K-dependent fine-tuning coefficient. The detailed fitting is shown in Appendix B. With 
them we can calculate F for each point and plot them in Fig. 3, from which we obtain: 

• A large A can lift the Higgs boson mass, but count against naturalness. In spite of a 
very small change in the numerical value of A, the available maximum my, shifts a few 
of GeVs. However, as A increases, the strong A— RGE effects become more significant 
and thus As is larger. If we slightly increase A to a value 0.65, my, can reach ~130 
GeV, but As typically lies above 100. Hence A being properly large instead of as large 
as possible, is favored by naturalness. 

• The new fine-tuning As tends to dominate over F m 2 , except in the large k region 

H u 

where the m 2 s self- reducing is significant. Nevertheless, a large k strongly prefers the 
pulling region rather than the pushing region as explicitly explained in Section II B 3. 
It also interprets the increasing of F with my,. 

• We are not going to present the fine-tuning in the pulling scenario in this work. It 
typically gives fine-tuning dominated by F m 2 > 60 which is worse than the pushing 
scenario. But As is never a problem as explained in the above. In this scenario, our 
results of the minimal fine-tuning > 60 is in agreement with the one found in Ref. [3] 
which uses the leading logarithm approximation. 

In short, we clarify the actual naturalness problem associated with the relatively heavy Higgs 
boson within the NMSSM (with a high mediation scale): We need a large A and a relatively 
small k to push the Higgs mass, but they incur a large fine-tuning in stabilizing the vacuum. 
Although \A\\ < \A K \ may also destabilize the vacuum, we have shown previously that it 
can be avoided at the price of a small fine-tuning. Our results generalize the conclusion of 
Ref. [31], where such kind of new fine-tuning is presented only in the PQ-limit. 

Now we turn to the Higgs boson collider phenomenology, focusing on the h — > 77 signal. 
At the ATLAS and CMS experiments, di-photon is the main excess for the hints of the SM- 
like Higgs boson with mass in the range 123-127 GeV. Interestingly, to some extent in the 
natural NMSSM the excess can be regarded as a prediction from the pushing scenario. To 
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FIG. 4: The ratio of h — > 27 between the NMSSM and SM predictions. Top figures: Pushing 
scenarios with A = 0.60 (left) and A = 0.62 (right). Bottom figures: Pulling scenarios with 
mg = 350 GeV and A t = -580 GeV. 

effectively push the Higgs boson mass, the singlet usually snatches a portion of component 
from H®. Furthermore, the required small tan (3 may reduce the bottom Yukawa coupling 
compared to the large tan /3 case. As a result, the h — > 27 signal is likely to be enhanced [32] 
mainly by reducing the branching ratio BR(/i — > 66). By contrast, in the MSSM to have 
such a heavy Higgs boson is already a rather tough task, and the di-photon excess is even 
more difficult to get [33]. Numerically, we adopt the method used in Ref. [15] to calculate 
the ratio of di-photon events in the NMSSM to the corresponding value in the SM: 

Rh( } = T(gg -> h) x BR(h -> 77) m 

From the NMSSMtools we can extract the decay widths ratio BR(/i — > jrffBR^hsM — > 77), 
and the Higgs production rates ratio T(gg — >■ h)/T(gg — > /ism), which is approximated by 
the square of the reduced couplings ratio Cg g = 9hggl 9h SM gg- We plot R h {^) in Fig. 4, from 
which some observations are made: 
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• The pushing scenario shows significant di-photon excess, especially the ATLAS and 
CMS data can be well fitted in the A = 0.62 case. From Fig. 4 it is reasonable to 
conclude: Naturalness prefers the Higgs boson with a smaller mass nih ~ 123 GeV 
and R h {^ r y) < 2, which favors the CMS results. In natural SUSY the light stop sector 
affects the di-photon rate from two aspects. On the one hand, it increases the Higgs 
production cross section via gluon fusion to a mount [14] 



To get the estimation we have used Eq. (39) and set tan — 2, \i — 200 GeV. On 
the other hand, stops decrease the width of Higgs decay to di-photon. The degree of 
decrease is about 20%, if the reduced coupling ratio between Higgs and vector bosons 
is quite close to 1. So the light stop sector only mildly increase the di-photon rate by 
~ 10%. Generically, the excess from light stops is about a few times 10%'s. But it 
is still very important to isolate this contribution, and we shall leave it for our future 
work. 

• By contrast, it is harder to achieve a relatively heavy Higgs boson and obvious di- 
photon excess in the pulling scenario. In fact, the excess tends to be slightly below 
the SM prediction. 

In the small A scenario, rrih is expected to be pushed totally by the mixing effect. Although 
here vacuum stability does not recur fine-tuning, we suffer a severer dependence of m/, on 
the stop sector. Consequently, the fine-tuning is expected to be worse. Additionally, the 
H® component in h is considerably reduced in most discussed in Section II B 4. 

So the diphoton excess is not so spectacular as in the large A scenario. Fig. 5 supports 
the conclusion: We have only scanned the regions near the most optimistic points, but the 
points showing significant di-photon excess are still rare. 

The recent hints favor the pushing scenario with a large A, but we are looking forward to 
the confirmative precise measurements of the Higgs boson properties to distinguish different 
scenarios. 

B. The neutralino LSP dark matter: a challenge from XENON100 experiment 

The existence of dark matter is an important evidence for new physics, maybe regarded 
as a triumph for SUSY which naturally provides the neutralino LSP as a WIMP dark matter 
candidate. So it is necessary to investigate the naturalness implication on the neutralino 
DM. By definition, the natural NMSSM (We focus on the most interesting case with a large 




(42) 
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FIG. 5: Di-photon rates versus the Higgs boson mass in the small A scenario. Left plot is for 
nig = 400 GeV and A t = 800 GeV while the right for m f = 450 GeV and A t = 1100 GeV. 

A) displays a light neutralino world: the small M1/2 and /x, and a light singlino s as well 
due to Ms = 2kv s < Xv s = \i in the NMSSM. Furthermore, a large A implies that the bino, 
Higgsino, and singlino are well mixed in the neutralino. Consequently, the large Higgsino 
component [43] makes the LSP strongly disfavored by the XENON100 experiment. The 
Higgs bosons Hi mediate the DM-proton interaction, which gives rise to the spin-independent 
(SI) cross section a p . At the zero-momentum transfer it is formulated as a p = A/x^f^/ir [35] 
with the reduced mass \i p = Mdm^p/ (Mdm + m p ) ~ m p , and 



fp=^(-^ L a) E \°^ a + O iHd tan/3 /« 
m Hi \vsmp J a ~j~^ 3 L 



« (0^0.144 + 0.206O^tan/3) ^-^J ^p, (43) 

where the first and second terms in the bracket respectively denote the up-type and down- 
type quark contributions. We have used f$ = 0.020 ± 0.004, = 0.014 ± 0.003, f$ = 
0.026 ± 0.005, f$ = 0.036 ± 0.008, and f^ n) = 0.118 ± 0.062 [36]. Because the up-type 
quarks give the dominant contributions, we get 

On the other hand, for a typical WIMP with mass ~ 20 - 100 GeV, the XENON100 
experiment upper bound is f p ~ 0.5 x 10" 8 GeV" 2 . So generically a p is about one order 
larger than the experimentally allowed bound, justified by Fig. 6. 

However, the LSP can still circumvent the exclusion in a subtle way. The point is that the 
coupling of the vertex #211-^2X1X1 can be suppressed considerably as a result of accidental 
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FIG. 6: Left: The LSP neutralino relic density under the WMAP constraint versus its mass. 
The annihilation rate is not problematic due to the large Higgsino component. Right: the LSP 
neutralino-proton spin-independent recoiling cross section under the XENON100 experimental con- 
straint versus its mass. It is strongly disfavored. Here A = 0.60. 

cancellation (We noticed that this point has also been mentioned in Ref [37].). Interestingly, 
in the pushing scenario the WMAP and XENON100 experiments together "predict" a LSP 
with mass in the vicinity of ~ 63 GeV. While in the pulling scenario the cancellation 

happens for the heavier neutralino, so it "predicted" a larger mass ~ 100 GeV, as shown in 
the right panel of Fig. 6. If we relax the naive WMAP bound via a non-standard thermal 
history of the Universe, the XENON100 experiment leaves more room for the neutralino 
LSP in the pushing scenario, as shown in Fig. 6. However, once the naturalness bound is 
further imposed, this case is again not favored. Fig. 7 shows that in the region allowed 
by the XENON100 experiment the degree of fine-tuning is typically below 1%, in both the 
pulling and pushing scenarios. Therefore, we may want to consider the other more natural 
LSP DM candidates, such as the sneutraino [38]. 

Some remarks are in orders. (A) For a small A, the neutralino LSP passing the constraints 
also possess a large Higgsino component. But this scenario accommodates no cancellation, 
and then hardly reconciles with the XENON100 experimental bound. (B) A few works [44] 
have studied the correlation between the Higgs and DM searchs as well [39]. However, in 
the NMSSM a relatively heavy Higgs boson, following naturalness criteria, has more direct 
impact on the neutralino LSP DM: It prefers the light neutralino world, which is however 
severely disfavoured by the XENON100 experiment. Of course, this correlation between the 
naturalness and LSP DM search can be removed, if we give up the hypothesis of gaugino 
mass unification. In that case the LSP can be bino-like and hence evade the constraint from 
the naturalness discussed here. (C) If we stick to gaugino mass unification, the scenario 
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with low mediation scale stands out again, where the LSP is the gravitino which obviously 



escapes the XENON100 experimental bound. 
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FIG. 7: The naturalness implications on the LSP neutralino. The arrow indicates the XENON100 
experimental bound. Here, we have fixed A = 0.60. 



IV. DISCUSSION AND CONCLUSION 

We considered the naturalness of mz and rrih in the MSSM and NMSSM via the bottom- 
up approach. For the GUT-scale SUSY breaking mediation with light stops, the gluino 
dominates the source of fine-tuning of mz, and the present LHC lower bound already gives 
Az < 2.5%. Also, the relatively heavy rrih near 125 GeV strongly favors the NMSSM in view 
of naturalness. Then we studied the Higgs boson mass for a few scenarios in the NMSSM: 
(1) A large A in the pushing scenario which has a heavy Higgs boson and an significant di- 
photon excess. However, this scenario has a new large fine tuning, owing to vacuum stability 
and large A— RGE effects. (2) A large A but in the pulling scenario, A ~ 0.7 and tan/3 ~ 2 
allows a heavy m^. The difference between this scenario and the above one is the absence 
of a lighter Higgs boson. And the di-photon excess may not be generated in this case. (3) 
A small A, the mixing effect is still able to give a heavy enough Higgs boson. However, it 
may not have the significant di-photon excess. Also, the naturalness status becomes worse 
in both the second and third cases. In all these scenarios the LSP neutralino DM is strongly 
disfavoured by the XENON100 experiment. Even if the LSP neutralino DM is fine, the 
naturalness still disfavors it. Note that the above analyses are based on the mSUGRA-like 
model with GUT-scale mediation, whereas for the low scale mediation the naturalness can 
be improved by about one order. 

There are still some open questions to realize the natural NMSSM. For example, how 
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light the stops that the LHC can tolerate is still an important problem in the light stop 
framework, although some preliminary works appeared [4]. We stress that although the 
discussion on the new fine-tuning associated with ms is confined to the conformal NMSSM 
in this work, it may applies to the general NMSSM where the cancellation discussion in this 
work is still required. 
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Appendix A: An approximated solution to the Higgs boson masses 

For further convenience, we present the neutral Higgs fields in a basis with explicit Gold- 
stone boson G° [40], 

H° u =v u + cos f3 + S 2 smf3) + A=(Pi cos/3 + G°sin/3), 

v 2 v 2 

H° d =v d + - r (-5'i sin/3 + S 2 cos/3) + -^(P x sin/3 - G° cos/3), 
V 2 v 2 

P 2 + iS 3 

S =v s + . (Al) 

Thus, the elements of the CP-even Higgs boson mass squared matrix Mj (in the basis 
(S 1 ,S 2 ,S 3 )) are given by 

(M|) n = M 2 + (m| - A V) sin 2 2/3, 

(M|) 12 = -i(m|-AV)sin4/3, 

(M|) 13 = -(Ml sin 2/3 + 2 A/™ 2 ) cos 2/3— , 

(M|) 22 = m| cos 2 2/3 + A V sin 2 2/3, 

(M|) 23 = -(4A\ 2 - Ml sin 2 2/3 - 2Act 2 sin 2/3) — , 
2 v s 

(M|) 33 = sin 2 2/3 + 4/€ 2 w 2 + kA k v s - ^Xkv 2 sin 2/3, (A2) 

where Ml = 2\v s (A\ + kv s ) / sin 2/3. The doublet block has been approximately diagonalized 
already, with (Mg)| 2 being the frequently quoted tree-level upper bound on the CP-even 
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Higgs boson mass square. For completeness, we also give the two CP-odd Higgs boson mass 
squared matrix elements: 



(M|) n = M 2 A , 

(M|) 12 = ^(M 2 sin2/3 - 6\kv 2 s )-, 



(M|) 22 = -(Mlsin2/3 + 6A^ s 2 )sin2^ 




(A3) 



Considering the 3x3 CP-even Higgs boson mass squared matrix structure, we can find 
the quite precise approximate eigenvalues analytically. First, we diagonalize the 13-block 
and get two eigenvalues: 



with the corresponding eigenvectors dominated by the first component. Now the heaviest 
eigenvalue is isolated, and the interesting two lighter eigenvalues can be obtained from the 
following effective 23-block: 



Thus, it can be solved analytically. 

Appendix B: The m| dependence of k 

m 2 s is sensitive to k in the large A limit. The fraction of the m| component in m|, 
denoted as C m |(K, A), is important for the fine-tuning calculations. For A = 0.62, we use 
the following function to fit C m 2 (k, A): 






(A4) 




(A5) 



C m 2 («, A = 0.62) « 0.57 + 0.54k - 5.42k 2 . 



(Bl) 



And for A = 0.60, the fitting function is 



C m 2 («, A = 0.60) « 0.62 + 0.28k - 4.43k 2 . 



(B2) 



Moreover, the fitting function plot is given in Fig. 8. 
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